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Abstract 

By numerically investigating the conservation law of the supercurrent, we con- 
firm the restoration of supersymmetry in Sugino's lattice formulation of the two- 
dimensional N = (2,2) supersymmetric SU(2) Yang-Mills theory with a scalar 
mass term. Subtlety in the case without the scalar mass term, that appears to ruin 
perturbative power counting, is also pointed out. 
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1 Introduction 



Nonperturbative study of supersymmetric gauge theory is of great general 
interest but in the context of lattice formulation, no compelling evidence of 
supersymmetry has been observed so far. Spacetime lattice is generally irrec- 
oncilable with supersymmetry and one must fine-tune coefficients of relevant 
and marginal operators so that supersymmetric Ward-Takahashi (WT) iden- 
tities are restored in the continuum limit. (An important exception is the 
four- dimensional M = 1 supersymmetric Yang-Mills theory (SYM) pQ; see 
Ref. [2] for a lot of effort went into numerical study of this system; see Ref. [3] 
for a recent attempt.) Recently, for two- and three-dimensional extended su- 
persymmetric gauge theories, lattice formulations that require no (or a little) 
fine tuning have been proposed [llSISlTISigiTOITTITSITglT^ITSITmrTITglFI It can 
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be argued that exact fermionic symmetries in these lattice formulations, com- 
bined with other lattice symmetries, (almost) prohibit relevant and marginal 
supersymmetry breaking operators to appear [ofT] . Supersymmetry is then 
restored in the continuum limit without (or with a little) fine tuning^] 

The aim of the present study is to test this scenario of supersymmetry restora- 
tion, in Sugino's lattice formulation of the two-dimensional M = (2, 2) SU(2) 
SYM [7115] . By a Monte Carlo simulation, we study a supersymmetric WT 
identity in the form of the supercurrent conservation law. For the reason elu- 
cidated in Sec. HI we introduce a supersymmetry breaking scalar mass term 
to the lattice formulation. The expected WT identity hence takes the form of 
a partially conserved supercurrent ("PCSC") relation in which the breaking 
term is proportional to the square of the scalar mass. We numerically con- 
firm the restoration of this PCSC relation in the continuum limit. Our result 
strongly indicates that the proposed scenario for the supersymmetry restora- 
tion is in fact valid and the target continuum theory (i.e., the two-dimensional 
M = (2, 2) SU(2) SYM with a soft supersymmetry breaking mass term) is re- 
alized in the continuum limit of the present lattice model. 



2 Preliminaries in the continuum theory 



The euclidean continuum action of the two-dimensional Af = (2, 2) SU (N c ) 
SYM, that is obtained by dimensionally reducing the Af = 1 SYM from four 
to two dimensions, is given b)0D 

S = - 2 J d 2 x tr \^F MN F MN + m T CT M D M m + # 2 } , (3) 



2 See also Refs. [26 27J for alternative approaches that do not use exact fermionic 
symmetries. 

3 In this paper, we adopt the representation 
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that corresponds to the conventional basis in the topological field theory |28|29] . 
See also Ref. 1301. 
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where Roman indices M and N run over 0, 1, 2 and 3, while Greek indices fi 
and v below run over only and 1. F MN are field strengths F MN = Dm An — 
dxA M + i[A M , A N ] and covariant derivatives D M are defined with respect to 
the adjoint representation, D M ^ = + z[Am, Here, in all expressions, 
it is understood that 82 — > and ^3 — * (dimensional reduction). We also 
define complex scalar fields, <fi = A2 + iAg and <fi = A2 — 1A3. The auxiliary 
field H is related to the auxiliary field H in Refs. [7f8] by H = H — z\Foi- We 
also introduce a soft supersymmetry breaking term 

1 r 

Smass = J d 2 X /i 2 tr { </></>} . (4) 

This term is "soft" in the sense that it does not introduce new ultraviolet 
divergences compared with the supersymmetric theory S. 

Action ([3]) is invariant under the super transformation, 5Am = ie T CTM^, 
5^/ = ^FMN^M^N^ + iHT 5 e and 5H = —ie T CT 5 TMF>M^ , where the parame- 
ter e has four spinor components, e T = — (e^°\ e, e). We define components 
of the super transformation by 

6 = £ (0)q(0) + £ (i)g(D + £ q (5) 

Also, the total action 5* + S mass possesses following global symmetries: The 
U(1)a symmetry, 

^ — > exp {aT 2 r 3 } \1>, 0^exp{2za}0, — > exp {— 2ia} 0, (6) 
the U{l)v symmetry, 

^ -y exp{iar 5 }*, (7) 

and a global Z 2 symmetry, 

^ -> *r 2 ^, 0^-0, -> -0. (8) 

The last one is a remnant of the reflection symmetry with respect to the M = 2 
direction of the four- dimensional theory before dimensional reduction. 

A Noether current associated with the supersymmetry of 5 (the supercurrent) 
is given by 

s„ = -^CTmTnT^tIFmn^} = (J^\J£\J^J,) T ■ (9) 

These four spinor components correspond to fermionic transformations in Eq. (J5]), 
Q(°\ Q and Q, respectively. If the auxiliary field H, instead of H, is re- 
garded as an independent dynamical variable as is the case below, F 01 = — F 10 
in this expression must be understood as — iH. 
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We will consider a correlation function of the supercurrent and another fermionic 
operator. As a lowest- dimensional gauge invariant fermionic operator, we take 

U = ^r,(T 2 tr{A 2 ^} + T 3 tr{A^}) = (xf , X£\ X„ X,f . (10) 

This is gauge invariant because the scalar fields A 2 and A3, as well as the 
spinor transform as the adjoint representation under two-dimensional gauge 
transformations. We in particular consider a "diagonal part" in the product 
of supercurrent (Q and operator ( ITUl) . That is, we consider following four 
correlation functions 

(Jf(^f(0)), (^(^(O)), (j»(x)X„(0)) , (J,(x)X u (0)) . 

(11) 

It is useful to note that these four correlation functions are interchanged each 
other under transformations (J7|) with a = tt/2 and (JED- Since Eqs. (J7|) and (jHJ) 
are symmetries of the present system, four correlation functions in Eq. (ITTj) 
coincide in the continuum theory. On the other hand, the present lattice for- 
mulation is not invariant under either Eq. (J7|) or Eq. ([8]) (while Eq. ([6]) is 
manifestly realized) and a "degeneracy" of these four correlation functions 
can be used as a useful measure for how we are close to the continuum. 

With the presence of the supersymmetry breaking term .Smass (HI), the super- 
current is not conserved and, defining 

/ = -2C (r 2 ti{A 2 ^} + r 3 ti{A^}) ee (y<°>, yw y, y) T , (12) 



we obtain the supersymmetric WT identity (the PCSC relation) in spinor 
components 
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if the regularization preserves supersymmetry. We emphasize that these lo- 
cal relations hold irrespective of boundary conditions. One can derive these 
relations in the functional integral by employing a local change of variables 
that does not "touch" the boundary. See Appendix [A] Therefore, these hold 
in particular with the antiperiodic temporal boundary condition for fermionic 
variables, although this boundary condition explicitly breaks supersymmetry 
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(for example, the energy spectrum would not be supersymmetric with this 
boundary condition). 

We now argue that, if the argument for the supersymmetry restoration in Ref. [7] 
is valid, relations (fT3l) - ffT6l) must be reproduced for x 7^ in the continuum 
limit of the lattice model. 

According to the argument of Ref. [7J, the lattice action in Ref. [7J provides 
a regularization, that becomes supersymmetric in the continuum limit, for all 
1PI correlation function of elementary fields. In particular, possible ultraviolet 
divergent functions, tadpoles and the scalar self-energy at the one-loop level, 
take the form (in the continuum limit) expected in the continuum theory (i.e., 
tadpoles vanish and no self-energy correction at zero external momentum). In 
this sense, the present lattice regularization is supersymmetric. This almost 
implies that supersymmetric WT identities such as Eqs. ffT3l) - ffT6l) hold in the 
continuum limit. 

Eqs. f|T5]) - ffTB]) . however, contain composite operators and a definition of com- 
posite operators in a lattice formulation is to a large extent arbitrary. We 
thus must be sure that lattice artifacts, when combined with ultraviolet diver- 
gences arising from these operators, do not modify the WT identities. Fortu- 
nately, in the present two-dimensional super-renormalizable system, operators 
in Eqs. flU]), (TTOT) and (fT2j) themselves are ultraviolet finite (i.e., no operator 
renormalization/mixing is required; this significantly simplifies consideration 
of the supersymmetry WT identity compared with four dimensions) Q Only 
possible ultraviolet divergence in Eq. (fT6l) . for example, arises from a one- loop 
diagram obtained by connecting elementary fields in and X v . However, this 
ultraviolet divergence can readily be avoided by setting two points x and to 
be separated!!] In fact, by analysing this one-loop contribution, one has 

d» (Mx)X u (0)) - !L (Y(x)X u (0)} 

= -i5 2 {x) (QX u (0)) + ^(N 2 C - l)(c Q - l)dj\x) (17) 

for Eq. (fT6j) and similar corrections for Eqs. (fT3]) - (jT5l) (the constant cq here 
can differ for each of relations, as Cq(o), Cq(i) and c^, depending on the regu- 
larization), assuming that the regularization for 1PI functions of elementary 
fields preserves supersymmetry. The constant cq in Eq. (JTTj) depends on the 
definition (or regularization) of composite operators and X v , but this de- 

4 Only possible operator that the supercurrent can mix with (in the one-loop level) 
is CT^tr-f'Q/} and this identically vanishes for the gauge group SU(N C ). 

5 This is impossible for one-point WT identities studied in Refs. |31|32|33|34j ; these 
one-point functions are thus subject of lattice artifacts and cannot directly be used 
to observe the supersymmetry restoration. 
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pendence on regularization disappears for x 7^ 



In summary, if the argument for the supersymmetry restoration in our lattice 
model is valid, relations f|T3l) - f|T6l) with x ^ must be reproduced in the 
continuum limit PI 

One might wonder why we worry about the supersymmetry restoration in 
the present model; after all, supersymmetry is explicitly broken by the scalar 
mass term (and by the boundary condition that we will adopt below). It is 
very important, however, to distinguish three different possible sources for su- 
persymmetry breaking in the present lattice model; the scalar mass term, the 
boundary condition and the lattice regularization itself. Our point is that the 
observation of WT identities (fl3"]) - (fIS]l enables us to isolate the last source 
of supersymmetry breaking. What we are talking about is whether this su- 
persymmetry breaking owing to the lattice regularization disappears in the 
continuum limit or not. 



3 Results of Monte Carlo study 

For details of the lattice formulation in Refs. [718], we refer the reader to the 
original references and Refs. [33|35] . The point is that the fermionic transfor- 
mation Q in Eq. ([5]) is nilpotent (up to gauge transformations) and action (j3J) 
can be written as a Q-exact form, as topological field theory [28|29] . These 
nilpotent Q and Q-exact action, and thus the invariance of the action un- 
der Q, can be realized even in lattice gauge theory [7f8"] . This sort of exact 
fermionic symmetry is realized also in lattice formulations of the present sys- 
tem in Refs. jlf5flOM12|[T3"] . Other fermionic symmetries, Q^°\ and Q, from 
which Eqs. (fT3|) - (|T5|) follow, are not preserved in the lattice formulation and 
the question is whether the invariance under these is restored in the continuum 
limit. 

We consider a finite two-dimensional rectangular lattice 

A = [x e al? | < x < = 2L, < x x < L} , (18) 

where a is the lattice spacing and the physical size L is fixed to be v^/glZI 

6 Incidentally, the last term of Eq. (|17p is not a genuine anomaly of supersymmetry. 
In fact, there exists a possible definition of that leads to cq = 1; see footnote [TT1 

7 Our argumentation here is not quite rigorous. For more satisfactory treatment, 
one first derives a WT identity associated with an appropriate (would-be) super- 
symmetry transformation in lattice theory and then shows all lattice artifacts to 
Eqs. (|13p -(|16 p with i/O vanish in the continuum limit. 

8 Note that in two dimensions, the gauge coupling g has the mass dimension 1. 
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Except the result in Fig. [7J below, all the results were obtained with the an- 
tiperiodic temporal (= Xq) boundary condition for fermionic variables (the 
aq-direction is always periodic). We found that this boundary condition leads 
to faster approach to the continuum and less noisy signals, compared with 
the periodic boundary condition. Recall that WT identities (fT3]) - (fT6l) in the 
continuum theory must hold irrespective of the boundary condition] 

Corresponding to mass term we add a scalar mass term 

H 2 
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9 2 x&A 



J2tr{<P(x)<P(x)} (19) 



to the lattice action of Ref. [8][_^j For lattice transcription of operators in 
Eqs. (JUJ), ffTUl) and (fT21) . we adopt a simple prescription that a field variable 
at a point x is replaced by the corresponding lattice field variable at x. Co- 
variant derivatives for the scalar field are replaced by the forward covariant 
differences, D^(f>(x) — > {U(x, /i)0(x + ajl)U(x, /i) _1 — 0(x)}/a 2 , where U(x,fi) 
are link variables. As noted in the preceding section, a precise way of lattice 
transcription of these operators should be irrelevant to our present analysis. 
That is, all ambiguities associated with the operator definition are integrated 
in the constants eg in Eq. ffTTl) (and Cg(o), Cg(i) and Cg) in the continuum limit 
and they do not appear for x ^ O P^l 

One of us (I. K.) developed a Rational Hybrid Monte Carlo (RHMC) |36j sim- 
ulation code (on the basis of the multi-shift [37] CG solver |38j ) for the present 



9 In the present lattice formulation, the invariance under the fermionic transfor- 
mation Q is manifest. Then, one can derive an exact WT identity of the form (|16|) 
in lattice theory with an appropriate definition of composite operators (see foot- 
note [TT]). The boundary xq = = j3, to which the antiperiodic boundary condition 
refers, can freely be shifted by change of fermionic variables. This also illustrates 
that the boundary condition is irrelevant to local WT identities such as ^3|) - ^6|) . 

10 We adopt a convention that all lattice field variables are dimensionless; the con- 
tinuum dimensionful scalar field <j> is replaced by <j)(x)/a on the lattice. 

11 By employing lattice perturbation theory and the Reisz power counting theorem, 
one finds that the above definition of Jn and X v leads to eg = 1 — ir in Eq. ()1T[) . 
Since the present lattice formulation has an exact fermionic symmetry Q, there 
exists an alternative natural definition of that exactly fulfills the corresponding 
WT identity even with finite lattice spacings. This improved Noether current differs 
from the above definition by terms of 0(a) and these terms improve the quantum 
property of the current. In particular, in the continuum limit, one has eg = 1 
in Eq. (|17p . Although this improved current has a desired property, here we do not 
adopt this definition for several reasons. Firstly, our prime objective is to see the 
restoration of symmetries other than Q; there is no particularly superior definition of 
Noether currents associated with Q^°\ and Q in the present lattice formulation. 
Secondly, we want to illustrate the idea that a precise form of the operators is 
irrelevant for relations (Tl3]) - p6]) with x ^ 0. 
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lattice system with the gauge group 577(2), by utilizing available libraries and 
programs |39|40f 41]. Some details of this simulation code have already been 
reported in Ref. [42J. (See also Fig. 3 of Ref. [13] that illustrates how the effect 
of dynamical fermions is properly included with this code.) 

We summarize on parameters in the present simulation: The parameter e for 



the admissibility [8] was fixed to be 2.6[H| The degree of the rational approx- 



imation in RHMC was typically 20-30. The multi time step acceleration [H] 
was used in the molecular dynamics (see Appendix [B]) . The time interval of 
one trajectory was fixed to be 0.5 and time steps of leapfrog were chosen such 
that the acceptance was greater than 0.8. We stored configurations in every 
10 trajectories. The autocorrelation was then estimated by jackknife analy- 
sis with binning. We discarded first 10000 configurations for thermalization 
and used subsequent configurations at every 50 configurations. In this way, we 
prepared uncorrelated configurations listed in Table [TJ 

In the present lattice formulation, the pfaffian of the Dirac operator resulting 
from the integration over fermionic variables, that is real positive in the con- 
tinuum theory, is generally complex owing to lattice artifacts. For parameters 
in the present simulation with the antiperiodic boundary condition, we found 
that the absolute value of the complex phase of the Dirac determinant is < 0.8 
radians. We took into account the complex phase of pfaffian by reweighting 
the half of the phase of the determinant. See Refs. [33,35j . We observed that, 
however, this reweighting had practically no notable effect. Also, all quantities 
presented below, that are real in the continuum theory, are generally complex 
with finite lattice spacings. The imaginary part is however generally small (and 
consistent with zero) and we will present only the real part for simplicity. 

For the correlation functions used in the following analyses, we took an av- 
erage over the whole lattice points to increase the effective number of con- 
figurations. For example, (Jq(x)Xq (Q)) is computed by taking the aver- 
age of (^Tq (v)Xq (z)) over all y G A and z G A with y — z = xq mod (3 
and i/i — z\ — x\ mod L are kept fixed 
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12 Although a precise choice of the parameter < e < 2\/2 should be irrelevant for 
results with small lattice spacings, we admit that we have not carried out systematic 
study on this point. As a simple consistency check, we performed a small experiment 
that measures the expectation value of the action density of the pure Yang-Mills 
part (Eq. (4.27) of Ref. [35J) for different values of e, e = 1.0, 2.0, and 2.6. The 
scalar mass squared is /x /g 2 = 0.25, lattice size is 6 x 6, and the lattice spacing 
is ag = 0.2357. The number of uncorrelated configurations is 100 for each e. We 
observed that the admissibility bound ||1 — U(x, 0, 1)|| < e (see Ref. [35]) was never 
exceeded in the molecular dynamics and the expectation values for each e coincided 
within statistical errors. 

13 With the understanding that the sign of the correlation function is flipped when 
y — z = x — (30, the correlation function is translationally invariant even with the 
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Table 1 

Set of configurations we used. 



v 2 /g 2 


lattice size 


ag 


number of configurations 


set label 


0.04 


12 x 6 


0.2357 


800 


1(a) 


0.04 


16 x 8 


0.1768 


800 


1(b) 


0.04 


20 x 10 


0.1414 


800 


I (c) 


0.25 


12 x 6 


0.2357 


800 


II (a) 


0.25 


16 x 8 


0.1768 


800 


II (b) 


0.25 


20 x 10 


0.1414 


800 


II (c) 


0.49 


12 x 6 


0.2357 


800 


III (a) 


0.49 


16 x 8 


0.1768 


1800 


III (b) 


0.49 


20 x 10 


0.1414 


1800 


III (c) 


1.0 


12 x 6 


0.2357 


800 


IV (a) 


1.0 


16 x 8 


0.1768 


1800 


IV (b) 


1.0 


20 x 10 


0.1414 


1800 


IV (c) 


1.69 


12 x 6 


0.2357 


800 


V(a) 


1.69 


16 x 8 


0.1768 


1800 


V(b) 


1.69 


20 x 10 


0.1414 


1800 


V(c) 



As a typical example of correlation functions, we plot in Fig. [T] correlation 
functions ffTTj) obtained by employing set IV (a) of Table [1] as a function of x 
(xi = L/2). Four correlation functions are quite well-degenerated, except at 
points near x = and its periodic images. (As noted, at x = and its periodic 
images, the correlation functions suffer from lattice artifacts associated with 
the operator definition.) We thus have an indication that the lattice spac- 
ing ag = 0.2357 is already rather close to the continuum, at least for the 
correlation functions with n 2 /g 2 = 1.0. 

It is instructive to plot the left-hand side of WT identities f TT5|) (TT^T) . If su- 
persymmetry is restored, the left-hand side must vanish except at x = and 
its periodic images; recall Eq. (|TT|) . In Fig. [2j we showed the left-hand side 
of Eq. ( fTBl obtained by set IV (c) of Table [TJ Here, we used the symmetric 
difference, d^f(x) = {f(x + aft) — f(x — afi)}/(2a), as a lattice transcription 
of the total divergence, because we found that its use considerably reduces the 
statistical error. We see that the left-hand side is almost zero everywhere (the 
absolute value is less than 0.05 in the central portion) and it is sharply peaked 
at the origin x = and its periodic images. This result strongly indicates that 

antiperiodic boundary condition. 
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Fig. 1. Correlation functions (jlll) with ji = v = along the line x\ = L/2 ob- 
tained by set IV (a) of Tabled (ji°\x)X^ (0))/g 2 (+), (J^ (x)X^ ] (0)) / g 2 (x), 
(Jo(x)X (0))/ 5 2 (□), (J (x)X (0))/g 2 (■). 

our reasoning in Sec. [2] is correct and, at the same time, that supersymmetry 
is certainly restored in the continuum limit. Recall that Eq. ( fT3l) is associated 
with the fermionic symmetry Q^> that is not preserved in the present lattice 
model with finite lattice spacings. 

For systematic quantification of the supersymmetry restoration, we consider 
the following four ratios: 



(21) 



(22) 



(r(°)(x)x ( S 0) (o) 



Y(x)X (oy / 

ffl (Mx)X (Q)) 
(Y(x)X (0)) 



(23) 



According to PCSC relation (|T3l) - (|T6l) . these ratios must become fi 2 /g 2 in 
the continuum limit, except at x = and its periodic images. We plotted 
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20 



Fig. 2. A bird's eye view of the left-hand side of Eq. (|13j) with v = (measured in 
a unit of g 3 ) on a 20 x 10 lattice. This was obtained by set IV (c) of Table [U 

ratio fl20l for different lattice spacings in Fig. [3j The statistical error in the 
ratio was estimated by jackknife analysis. This plot is for n 2 /g 2 = 1.0 and 
in the continuum limit the points should lie on the dotted line. We see this 
expected tendency, while a deviation near xq = = (3 can be understood 
as the effect of approximate delta functions at x = and periodic images 
elucidated above. Also, again, we see that the lattice spacing ag = 0.2357 
is already rather close to the continuum for fj, 2 /g 2 = 1.0, because the WT 
identity is fairly restored with this lattice spacing. 

To quantify the restoration of supersymmetry, while eliminating the effect of 
lattice artifacts existing around x = and its periodic images, we adopted the 
following procedure. First, we take a cylindrical region C in A, C = {x G A | 
j3/2 — Ai < x < (3/2 + A 2 , < xi < L). We then apply a x 2 -^ by a constant 
to ratios (|20|) - (|23|) obtained for x G C. We change Ai and A 2 such that the 
number of points in the region C becomes maximum while keeping the \ 2 P er 
one degree of freedom (dof ) reasonably small. The results of this procedure are 
summarized in Table [2] and in Fig. HI In Table [21 columns indicated by "used" 
show the number of points contained in the region C determined as above. 
In Fig. |U the x 2 -fitt e d values of ratio (|20|) obtained by the above procedure 
are plotted. As shown in the figure, these values were then used for a linear 
^-extrapolation to the continuum. 

Figure is the main result of this paper. The continuum limit of ratios (12"U|) - 
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Fig. 3. Ratio (|20j) along the line x\ = L/2. These were obtained by set IV (a) (+), 
set IV (b) (x) and IV (c) (□), respectively. 

( ]23l) obtained by the above procedure is plotted as a function of the parame- 
ter \J? I g 2 . The result is consistent with the straight line, that is a prediction 
of the supersymmetric WT identities (the PCSC relation). 

Some care should be paid for the interpretation of Fig. [5j the plot shows only 
statistical errors. There might be rather large systematic errors associated 
with the fitting procedure especially for the smallest mass /i 2 /g 2 = 0.04 case 
(set I). Figure M shows ratio (|2"U|) for set I and we see that a flat region is 
narrower compared with Fig. [31 that is for n 2 /g 2 = 1.0. In fact, as Table [2] 
shows, the numbers of points we used in the fit are rather fewer for set I 
than those for other sets. From Fig. [6j however, the systematic errors for 
set I would be at most 0.05, thus the result is still consistent with the PCSC 
relation even if this systematic error is taken into account. For other values 
of the mass parameters fi 2 /g 2 , the behavior of ratios is more or less similar to 
that in Fig. [3] and the numbers of points we used in the fit appear sufficient; 
so we do not expect large systematic errors. 

Therefore, from the agreement with the theoretical expectation in Fig. we 
infer that supersymmetry, that is broken in the present lattice regularization 
with finite lattice spacings, is certainly restored in the continuum limit J 14 



As a simple check for that nontrivial loop corrections are really important in 
our results, we repeated above analyses in the quenched approximation. We con- 
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Table 2 

The quality of the x 2_n t °f ratios f|20|) — (|23|) by a constant. 



set label 


Eq. 


csi 


Eq. 


(EH) 


Eq. 




Eq. 


(BSD 




used 


X 2 /dof 


used 


X 2 /dof 


used 


X 2 /dof 


used 


X 2 /dof 


1(a) 


12 


1.533 


12 


0.7085 


18 


0.7947 


12 


1.062 


1(b) 


24 


0.9475 


48 


0.9423 


48 


0.6633 


24 


0.5073 


1(c) 


20 


1.639 


20 


1.158 


70 


1.413 


20 


1.486 


II (a) 


18 


0.9632 


18 


0.4774 


18 


0.7021 


18 


0.8510 


II (b) 


16 


0.8439 


24 


0.9131 


16 


0.6593 


16 


0.8143 


II (c) 


50 


0.9912 


50 


1.047 


60 


0.9219 


50 


0.9078 


III (a) 


12 


0.5984 


18 


0.6636 


18 


0.9439 


18 


0.5896 


III (b) 


16 


1.676 


32 


1.347 


16 


1.441 


40 


1.416 


III (c) 


80 


0.9306 


80 


0.9628 


70 


0.9690 


80 


0.7924 


IV (a) 


12 


0.9697 


18 


0.8239 


12 


1.053 


12 


1.043 


IV (b) 


24 


0.9034 


32 


0.9519 


24 


0.8340 


24 


0.7721 


IV (c) 


50 


0.9871 


40 


0.9419 


50 


0.9792 


40 


0.9535 


V(a) 


18 


0.6278 


18 


0.9550 


12 


0.6992 


18 


0.6746 


V(b) 


32 


0.9747 


32 


0.9105 


16 


0.7900 


32 


0.6894 


V(c) 


60 


1.469 


80 


1.326 


70 


1.423 


70 


1.346 



The antiperiodic boundary condition explicitly breaks supersymmetry while 
the periodic one preserves supersymmetry, so it is certainly of interest to carry 
out simulations (for the energy spectrum, for example) with the latter bound- 
ary condition. A typical example of correlation functions with the periodic 
boundary condition is depicted in Fig. [71 The degeneracy of four correlation 
functions is not quite realized. This indicates that the spacing ag = 0.2357 
is not yet so close to the continuum for fi 2 /g 2 = 1.0 when we use the peri- 
odic boundary condition. In addition to this, we found that signals with the 
periodic boundary condition are generally rather noisy and, when translated 
to the ratio like Fig. [3], the statistical errors are too large for a reliable fit. 
For these reasons, we postpone a detailed study of cases with the periodic 
boundary condition to a future work. 



firmed that in fact the results significantly differ from those obtained above and the 
expected super symmetric WT identities are not restored. 
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0.05 



0.1 



0.15 



0.2 



0.25 



ag 



Fig. 4. The x 2 -fitted values of ratio (f20l) as a function of the lattice spacing ag. 
Linear ^-extrapolations to the continuum a = are also shown. Set I (+), set 
II (x), set III (□), set IV (■) and set V (())■ 

4 Subtlety in the supersymmetric /i 2 = case 



In the above analyses, we have introduced the scalar mass term (II 9 1) that 
generally suppresses the amplitude of scalar fields. One is of course interested 
in the original two-dimensional M = (2, 2) SYM that does not contain such a 
supersymmetry breaking term. We now explain why we had to introduce the 
scalar mass term in our simulation. 

Figure M shows correlation functions ( {TTi) without the scalar mass term. The 
lattice spacing is ag = 0.2357. One sees that the degeneracy among four 
correlation functions is enormously violated and thus it appears that we are 
quite far from the continuum theory. Even if we decrease the spacing to ag = 
0.1768 (Fig. [2]), the crude feature does not change and the degeneracy is not 
restored. This behavior prompts us to draw a conclusion that the degeneracy 
is not restored even in the continuum limit. 

Moreover, from a closer look at Figs. [S] and EH it appears that the non- 
degeneracy cannot be understood as a breaking of a certain symmetry in the 
continuum theory. In these figures, only (J'q°\x)Xq°\o)) (indicated by +) is 
quite different from others while other three are almost degenerated. As noted 
around Eq. (II ip . these four correlation functions in the continuum theory are 
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2.5 



1.5 - 



1 - 



0.5 




0.5 



1.5 



Fig. 5. The continuum limit of x 2 -fitted values of ratios (|20j) - (l23p as a function 
of /J, 2 /g 2 . Eq. £11 (+), Eq. (JZI]) (x), Eq. ([22]) (□) and Eq. £23]) (■). The straight 
line is a prediction of the PCSC relation. The plot for fJ, 2 /g 2 < 0.3 is magnified in 
the small window. 



related each other by certain global symmetriesJ 15 1 For example, (j7"o(x)X (0)) 
and (j7o(a;)Xo(0)) are related by transformation (J7J) with a = tt/2. The fact 
that these two functions are quite degenerated in the figures suggests that 
this global discrete symmetry is well restored with the present lattice spac- 
ings. However, (J'q°\x)Xq°\o)} and (J'q 1 \x)Xq 1 \o)) are not degenerated at 
all, although these two are also related by this symmetry in the continuum 
theory. Something quite strange seems happening. 



Another observation we made is that if we modify a term in the lattice Dirac 
propagator that originates from the Yukawa interaction, by an 0(a) amount 
(in the process of measurement), the effect of the modification is tremendous 
and four correlation functions become almost degenerated; the effect appears 
to be 0(1). 

We suspect that the above strange behavior in the absence of the scalar 
mass term is caused by very large expectation value of scalar fields along flat 
directions — continuous set of minima of the classical scalar potential. Here, 
by "very large", we mean the lattice cutoff scale, 0(1/ a). In fact, as Fig. 3 
of Ref. [12] shows, the (gauge invariant) amplitude of scalar fields can be as 



Note that the boundary condition does not affect these global bosonic symmetries. 
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Fig. 6. Ratio (|20l) along the line xi = L/2. These were obtained by set I (a) (+), 
set I (b) (x) and I (c) (□), respectively. 

large as ~ 40/ a, with the antiperiodic boundary condition] 16 

Such a very large expectation value could ruin perturbative power counting, 
in which the operator a<p, for example, is regarded as 0(a). That is, if the 
expectation value of the scalar field (f> is as large as 0(1 /a), the combina- 
tion a<j) would have to be regarded as 0(1). This phenomenon affects also 
the argument for the supersymmetry restoration. For example, the opera- 
tor Q(ati{(f)ip IM }) = atr{7yi/v} + a tr {(j)iD^(f)} is invariant under the gauge, 
U(1)a and Q transformations — manifest symmetries of the present lattice 
formulation — and thus could radiatively be induced (in the one-loop level). 
This is an irrelevant operator in the usual sense but this could be "marginal" 



On the other hand, we numerically observed that the amplitude does not grow 
so much (< 1/a) with the periodic boundary condition. This seems strange at first 
glance because the periodic boundary condition does not break supersymmetry and 
one may expect that the flat directions are not lifted upon radiative corrections 
when supersymmetry is preserved. This phenomenon could be understood if we 
consider an effective potential for scalar zero modes that is obtained by integrating 
out all other modes. With the antiperiodic boundary condition, there is no fermion 
zero mode and the degrees of freedom that integrated out are balancing between 
bosons and fermions. On the other hand, with the periodic boundary condition, 
fermionic zero modes are left unbalanced and they can induce nontrivial corrections. 
See Refs. |45|46j for related observations. 
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Fig. 7. Correlation functions (lllj) with v = along the line x\ = L/2. The lattice 
size is 12 x 6 and ag = 0.2357. fi 2 j g 2 = 1.0. The boundary condition is periodic. The 
number of configurations is 800. (J^> \x)xf ] '(0)) / ' g 2 (+), (J^ (x)X^ (0)) /g 2 (x), 
(J (x)X (0))/g 2 (□), (J o (x)X (0))/g 2 (■)• 



if the scalar field <fi is regarded as 0(1/ a). This 0(1) operator, moreover, is 
not invariant under Q^°\ and Q transformations, thus invalidates the 
argument for the supersymmetry restorationP^I 



It is thus conceivable that the target theory, the two-dimensional M = (2,2) 
SU(2) SYM, is not realized in the continuum limit of the present lattice for- 
mulation, unless we supplement the scalar mass term (or something that sup- 
presses the amplitude of scalar fields). This point is further discussed in the 
next section. Also, in light of this observation, a study on the dynamical 
supersymmetry breaking in the two-dimensional M = (2, 2) SU (2) SYM in 
Refs. P51H7] must be reconsidered [48J. 



7 One might think that the above operator is prohibited to appear in the effective 
action owing to the lattice rotation and reflection symmetries. We should note, 
however, that the present lattice formulation does not possess such symmetries. 
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Fig. 8. Correlation functions (jlip with /j, = v = along the line x\ = L/2. 
/j?/g 2 = 0. The lattice size is 12 x 6 and ag = 0.2357. The boundary condi- 
tion is antiperiodic. The number of configurations is 800. (J'q°\x)Xq \o)} / g 2 (+), 
(jU( x ) X W(0))/g 2 (x), (J (x)X (0))/g 2 (□), (Jo(x)X (0))/9 2 (■)• 

5 Conclusion and discussion 



In this paper, for the first time to our knowledge, the restoration of supersym- 
metry in a lattice formulation of a supersymmetric gauge theory was clearly 
observed. The PCSC relation in Fig. O first and foremost, can be taken as 
a solid basis for the lattice model in Ref. [8] to be used for evaluation of 
various quantities in the two-dimensional M = (2, 2) SYM with a supersym- 
metry breaking scalar mass term. Also, our result for a two-dimensional model 
demonstrates validity of general reasoning for the supersymmetry restoration 
on the basis of lattice symmetries and power counting. It is quite conceivable 
that, therefore, various lattice formulations of supersymmetric gauge theory 
that are based on similar reasoning work as they are aimed at. 

The above statement is for the lattice model in which there is no flat direction 
in the classical potential. The lattice model for the original M = (2, 2) SYM 
without the scalar mass term [8] possesses flat directions and, as we have 
argued in Sec HJ it is quite plausible that the target theory is not obtained by 
the continuum limit. Thus, we close this paper by summarizing the situation 
concerning the M = (2, 2) SYM without supersymmetry breaking terms: 



18 



1 

0.8 
0.6 
0.4 
0.2 



-0.2 
-0.4 
-0.6 



v X ^ X x ^ 

*i ■ ■ ■ ■ ■ x 



-m ~ X 



1 



+ + 



+ + 



+ + + 



+ + 



+1 



0.5 



1.5 

x g 



2.5 



Fig. 9. Correlation functions (jlip with = u 
fj?/g 2 = 0. The lattice size is 16 x 8 and ag 



along the line x\ = L/2. 
0.1768. The boundary condi- 



tion is antiperiodic. The number of configurations is 800. (J'q \x)Xq V> (0)) / g 2 (+ 



(J^(x)X^(0))/g 2 (x), (J (x)X (0))/g 2 (□), (J (x)X (0))/g 2 (■). 

(1) One may first introduce the scalar mass term and then take the limit [i 2 — > 
0. Physically, as discussed in Ref. Ej. if the mass \i is sufficiently small com- 
pared with the inverse of the system size, l/L, the effect of the breaking mass 
would be practically negligible because one cannot observe the wavelength 
longer than the size of the "universe" . This provides a possible way of defining 
the M = (2, 2) SYM. This route of definition would be mandatory in lattice 
formulations in Refs. [51ll2f24j . in which a supersymmetry breaking scalar mass 
term must be introduced to stabilize the lattice spacing. Our result in FigO 
suggests that this prescription works in the formulation of Ref. [8], because we 
have observed the restoration of supersymmetry even for /i 2 < l/L 2 = 0.5g 2 . 

(2) One may carry out Monte Carlo simulations without introducing the su- 
persymmetry breaking mass term. This is possible for example in the formula- 
tions of Refs. [5f2l)] . In this approach, however, there will be a subtle problem 
we encountered in Sec. H] that O(a) lattice artifacts seem to be amplified to 
0(1) by 0(1/ a) expectation value of scalar fields along flat directionsJ_|£j Our 
observation suggests that it is generally quite difficult to realize supersym- 



18 In the large N c limit, there is a possibility that the expectation value of scalar 
fields along flat directions is suppressed and one can evade the problem without 
breaking supersymmetry. See Ref. 
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metric theories with flat directions as a continuum limit of a lattice model, 
without suppressing the amplitude of scalar fields, because lattice formulation 
generally cannot be free from 0(a) discretization errors. 

A profound question is, however, whether above prescription (1) really pro- 
vides a "correct" definition of the target two-dimensional theory or not; the 
existence of the vacuum moduli is unlikely in two dimensions while the pre- 
scription would enforce scalar fields to localize around the origin. Unfortu- 
nately, we do not have the convincing answer to this question at present. 

We would like to thank Masanori Hanada, Hikaru Kawai, Martin Liischer, 
Hideo Matsufuru and Fumihiko Sugino for enlightening discussions. Our nu- 
merical results were obtained using the RIKEN Super Combined Cluster 
(RSCC). I. K. is supported by the Special Postdoctoral Researchers Program 
at RIKEN. The work of H. S. is supported in part by a Grant-in-Aid for 
Scientific Research, 18540305. 



A Derivation of supersymmetric WT identities ( 1131) — ( JT6I) 

In this appendix, supersymmetric WT identities ( fT3l) -( fT6l) are derived with 
emphasize on their independence of boundary conditions. We start with the 
expectation value of operator ffTO]) : 



where Z is the partition function and d\i denotes a measure of the functional 
integral. We then consider a certain infinitesimal variation of integration vari- 
ables 5', A M -> A M + S'A M , ^ -> ^ + 5'^ and H -> H + 5'H. Since the 
functional integral itself is independent of any relabeling of integration vari- 
ables, we have the identity for any variation 5', 



provided that the measure d[i is invariant under the variation 5'. 

Now, as a particular form of 5', we take the super transformation with a re- 
placement e — > e(x), where e(x) is a Grassmann-odd spinor function with a 
finite support that does not overlap with the boundary; 5'A M = ie 1 (x)CT M ^ ', 



5'V = ^FMNT M T N e(x)+iHT 5 e(x) and 5'H = -ie T (x)CT 5 T M D M ^+e T (x)CT D 1 ^- 






(A.2) 
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e 1 \x)CY]_Dq^ . Under this variation 



5' (S + S mass ) =i J d 2 x | ~d^ T (x)s )1 (x) - ^e T (x)f(x) 

-^TaTbT^Fab + 2F A F tlA + r M r 5 #j *| , 

(A.3) 

where combinations s M and / are given by Eqs. (jHJ) and ( TT2T) . respectively. In 
this expression, the second line identically vanishes for any boundary condition 
because e(x) = at the boundary. Similarly, in the first line, we may perform 
integration by parts neglecting boundary terms again because e(x) = at the 
boundary. Finally, setting e(x) to be the delta function (times a Grassmann- 
odd constant spinor), we have Eqs. (JX3J) — (JXSD as particular cases of Eq. (1A.2j) . 
The assumption that the measure dfi is invariant under 5' corresponds to, in 
the present formal treatment, the assumption that regular izat ion preserves 
supersymmetry. 



+ \j d 2 xd^iL T {x)c 



B Multi time step acceleration in our simulation 



Let Sb be the action consists only of gauge and scalar fields and SW the 
action of pseudofermions. SW is bi-linear in pseudofermions. In the molecular 
dynamics, if there exists a definite hierarchy between force originates from Sb 
and that from Spp, and if the latter is smaller than the former, one may reduce 
computational cost for the latter (it is expensive requiring inversion of a Dirac 
operator) by using different time steps for each of force. This is the multi 
time step acceleration [H]. In this scheme, the time-evolution operator in one 
trajectory At = ne is symbolically written as 

T(S PF , 6/2) {[T(S B , e/k)] k T(S PF , e)}"^ [T(S B , e/k)f T(S PF , e/2), (B.l) 

where T(S, e) denotes the time-evolution operator with respect to the action S 
in the time step e. That is, one uses a fc-times larger time step for force 
from Spp- 

In our Monte Carlo simulation, as shown in Fig. IB. 11 we found that typically 
force from Sp F is several times smaller than that from Sb and the scheme 
is in fact very efficient. In the example in the figure, our choice was k = 3 
and n = 6 in Eq. (IB. II) (we set At = ne = 0.5). Since the variation of the 
force is large (typically the same order of magnitude as the average itself), we 
chose a smaller value of k than naively expected from the figure. As a general 
tendency, when a lattice spacing becomes smaller, force from Sb becomes 
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Fig. B.l. Force for gauge and scalar fields (Fb) and that for pseudofermions (Fpp), 
for a 12 x 6 lattice with ag = 0.2357 and n 2 /g 2 = 0.25. The values are averaged over 
each trajectory and plotted every 100 trajectories. Force from 5b is about 6 times 
larger than that from <Spf- 



larger and, correspondingly, we could use larger k. In fact, we used k 
n = 6 for ag = 0.1768, and k = 8, n = 6 for ag = 0.1414. 



5. 



References 



[1] G. Curci and G. Veneziano, Nucl. Phys. B 292 (1987) 555. 

[2] I. Montvay, Int. J. Mod. Phys. A 17 (2002) 2377 [arXiv:hep-lat /0112007] . 

[3] J. Giedt, R. Brower, S. Catterall, G. T. Fleming and P. Vranas. larXiv:0807.2032l 
[hepdat] . 

[4] D . B. Kaplan, E . Katz and M. Unsal, JHEP 0305 (2003) 037 
[arXiv:hepdat /020601~9] . 

[5] A. G. Cohen, D. B. Ka plan, E. Katz and M. Unsal, JHEP 0308 (2003) 024 
[arXiv:hepdat /0302017] . 

[6] A. G. Cohen, D. B. Ka plan, E. Katz and M. Unsal, JHEP 0312 (2003) 031 
[arXiv:hepdat/0307012|. 



[7] F. Sugino, JHEP 0401 (2004) 015 arXiv:hepdat/0311021 1. 



22 



[8] F. Sugino, JHEP 0403 (2004) 067 [arXiv:hep-lat/0401017] . 
[9] F. Sugino, JHEP 0501 (2005) 016 [arXiv:hep-lat/0410035] . 



[io; 

[11 

[12 

[13! 
[14 

[15 

[16 
[17! 
[18! 
[19 
[20! 

[21 

[22! 

[23! 

[24 
[25 

[26! 
[27! 

[28! 
[29! 
[30 



S. Catterall, JHEP 0411 (2004) 006 [arXiv:hep-lat /0410052| . 

D. B. Kaplan and M. Unsal, JHEP 0509 (2005) 042 [arXiv:hep-lat /0503039| . 

A. D'Adda, I. Kanamori, N. Kawamoto and K. Nagata, Phys. Lett. B 633 
(2006) 645 |arXiv:hep-lat/0507029| . 

F. Sugino, Phys. Lett. B 635 (2006) 218 |arXiv:hep-lat/0601024] . 

M. G. Endres an d D. B. Kaplan, JHEP 0610 (2006) 076 
[arXiv:hep-lat /0604012] . 

A. D'Adda, I. Kanamori, N. Kawamoto and K. Nagata, Nucl. Phys. B 798 
(2008) 168 [arXiv:0707.3533l [hep-lat]]. 

K. Nagata and Y. S. Wu, arXiv:0803.4 339 [hep-lat]. 

S. Matsuura, arXiv: 0805 .44911 [hep-th], 

F. Surino. larXiv:0807.2683 l [hep-lat]. 



M. Unsal, JHEP 0610 (2006) 089 |arXiv:hep-th/0603046|. 

P. H. Damgaard and S. Matsuura, JHEP 0707, 051 (2007) |arXiv:0704.2696l 
[hep-lat]]. 

T. Takimi, JHEP 0707 (2007) 010 [arXiv: 0705 .38311 [hep-lat]]. 

P. H. Damgaard and S. Matsuura, JHEP 0708, 087 (2007) }arXiv:0706.3007l 
[hep-lat]]. 

P. H. Damgaard and S. Matsuura, JHEP 0709, 097 (2007) jarXiv:0708.4"T29l 
[hep-lat]]. 

S. Catterall, JHEP 0801 (2008) 048 [arXiv:0712.2532l [hep-th]]. 

P. H. Damgaard and S. Matsuura, Phys. Lett. B 661, 52 (2008) 
[arXiv:0801. 2936 [hep-th]]. 

H. Suzuki and Y. Taniguchi, JHEP 0510 (2005) 082 [arXiv:hep-lat /0507019] . 

J. W. Elliott and G. D. Moore, PoS LAT2005 (2006) 245 [JHEP 0511 (2005) 
010] |arXiv:hep-lat/0509032| . 

E. Witten, Commun. Math. Phys. 117 (1988) 353. 
E. Witten, Int. J. Mod. Phys. A 6 (1991) 2775. 

J. Kato, N. Kawamoto and Y. Uchida, Int. J. Mod. Phys. A 19 (2004) 2149 



|arXiv:hep -th/0310242|. 
[31] S. Catterall, JHEP 0603 (2006) 032 |arXiv:hep-lat70 602004 1 . 



23 



[32] S. Catterall, JHEP 0704 (2007) 015 [arXiv:hep-lat/0612 008|. 

[33] H. Suzuki, JHEP 0709 (2007) 052 |arXiv: 0706 .13921 [hep-lat]]. 

[34] H. Fukaya, I. Kanamori, H. Suzuki and T. Takimi, PoS LAT2007 (2007) 264 
[arXiv:0709.4076l [hep-lat]]. 

[35] I. Kanamori, F. Sugino and H. Suzuki, Prog. Theor. Phys. 119 (2008) 797 
[arXiv:071 1.21321 [hep-lat]]. 

[36] M. A. Clark, A. D. Kennedy and Z. Sroczynski, Nucl. Phys. Proc. Suppl. 140 
(2005) 835 l arXiv:hep-lat/0409133| . 

[37] A. Frommer, B. Nockel, S. Giisken, T. Lippert and K. Schilling, Int. J. Mod. 
Phys. C 6, 627 (1995) [arXiv:hep-lat/9504020] . 



[38] B. Jegerlehner, arXiv:hep-lat/9612014 



[39] M. Di Pierro, Comput. Phys. Commun. 141 (2001) 98 [arXiv:hep-lat /0004007] . 



[40] M. Di Pierro and J. M. Flynn, PoS LAT2005 (2006) 104 
[arXiv:hep-lat /0509058 1 . 



[41] M. A. Clark and A. D. Kennedy, http://www.ph. ed.ac.uk/~ mike/remez, 2005 



[42] I. Kanamori. TarXiv: 0809 .06551 [hep-lat]. 
[43] I. Kanamori, arXiv: 0809 .06461 [hep-lat]. 

[44] J. C. Sexton and D. H. Weingarten, Nucl. Phys. B 380 (1992) 665. 

[45] H. Aoki, S. Iso, H. Kawai, Y. Kitazawa and T. Tada, Prog. Theor. Phys. 99 
(1998) 713 [arXiv:hep-th/9802085| . 

[46] S. Catterall and T. Wiseman, JHEP 0712 (2007) 104 [arXiv:0706.35T8l [hep- 
lat]]. 

[47] I. Kanamori, H. Suzuki and F. Sugino, Phys. Rev. D 77 (2008) 091502 
[arXiv:071 1.20991 [hep-lat]]. 

[48] I. Kanamori, work in progress. 

[49] K. N. Anagnostopoulos, M. Hanada, J. Nishimura and S. Takeuchi, Phys. Rev. 
Lett. 100 (2008) 021601 |arXiv:0707.4454l [hep-th]]. 



24 



